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Abstract. To approximate solutions of a linear differential equation, we project, 
]~~] via trigonometric interpolation, its solution space onto a finite-dimensional space 

of trigonometric polynomials and construct a matrix representation of the dif- 
^jQ ferential operator associated with the equation. We compute the ranks of the 

matrix representations of a certain class of linear differential operators. Our 
numerical tests show high accuracy and fast convergence of the method applied 
to several boundary and eigenvalue problems. 
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r-i 1. Introduction 

^ In this paper, wc use trigonometric interpolation to approximate solutions of a 
differential equation Au — /, whose differential operator A with domain D{A) is 

> a formal polynomial of operators {l,x, and / e Range(^). A solution u is 

^ projected onto the space 7^ of L-periodic trigonometric polynomials of degree n. 

^ The projection T[u] , defined as a trigonometric interpolant of u, is identified with 

od 

Q a vector u of its vahies at = 2n + l partition points via an isomorphism ijj, where 

7—i n G N. The operator A is represented by an A'" x A'" matrix A defined implicitly 

by Av — ipT [ATlv]] for all v e D(A). The original equation is approximated by 

a system of linear equations Au — f , where f = i^Tlf]. 
d 
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In sect. El we prove that ii A = ao + ai j- + . . . + as4^ , then 



rank A = N + \ sign ao \ — 2m — 1 , 



where m > is the number of solutions, in the set {1,2, . . . ,n}, of a system of 
polynomial equations dependent on the set of coefficients {ajj^^o- 

Our method employs ideas similar to those of the method of discrete approx- 
imations of differential operators that uses Lagrange interpolation. The latter 
approach was developed and tested by Calogero O |6] to solve eigenvalue problems 
for ordinary and partial differential operators. Using this method, Mytropolski, 
Prykarpatsky, and Samoylenko created a general algebraic-projection framework 
and furthered the applications to Mathematical Physics |10j . 

Using Lagrange interpolation by algebraic polynomials, Calogero introduced a 
finite-dimensional matrix representation K G R^^^ of a differential operator /C 
that is a formal polynomial of |l,x, A differential equation )Cu = /, where 
/ G Range(/C), is then approximated by a system of linear equations Ku = f, 
where f G is a finite representation of /. Bihun and Lustyk estimated the 
approximation error of the method [31 E]. 

Bihun and Prytula proved that the rank of the representation K G M^^^ of a 
differential operator 



is N — k [2]. Therefore, completion of the discrete approximation Ku = f with 
initial or boundary conditions may lead to an inconsistent system of linear equa- 
tions. In [21 15], this obstacle was overcome by employing substitutions. Also, the 
authors recognized that Lagrange approximation is less effective for approximating 
periodic solutions of differential equations. Our study of a discrete approximation 
A G M^^^ of /C, obtained using trigonometric interpolation, aims to provide a 
class of differential operators whose representations' ranks are compatible with 
the dimension of the solution space of the differential equation }Cu = f. 



^ — ^k——f: + O-k+l 




+ • • • + Q-k+p 
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In two short communications [3 |8], Calogero explored the idea of approxima- 
tion of differential operators using trigonometric interpolation. He proposed two 
matrix representations of the differential operator ^ and calculated their eigen- 
values. It was proven that if a trigonometric polynomial u of degree n solves a 
differential equation Au = 0, where A = J2'j=o'^j'^^ then the solution of the 
approximate problem Au = can be used to recover the exact solution u, where 
A G M^^^ is a matrix representation of A constructed using trigonometric in- 



terpolation (c/. theorem [5] in subsection 5.2). To our knowledge, the method 
of discrete approximation of differential equations via trigonometric interpolation 
has not been implemented and tested. Our numerical tests show high accuracy 
and fast convergence of the method applied to several boundary and eigenvalue 
problems for differential equations. 

2. Trigonometric interpolation 

Definition 1. A trigonometric polynomial T{x) of period L and degree n is defined 
to be 



(1) T{x) = c + ^ttfc sin + cos ^^^^ 



k=l 

where c, ai, an, bi, bn are real numbers and a'^ + b"^ > 0. 

We will refer to trigonometric polynomials as simply polynomials. Should the 
need arise to discuss ordinary polynomials, they will be referred to as algebraic 
polynomials. We denote the space of all L-periodic trigonometric polynomials of 
degree n by T^. 

In [21 |5l [To] , Lagrange interpolation is used to construct finite-dimensional rep- 
resentations of operators |l,x, We construct similar representations using 
trigonometric polynomials to solve differential equations and eigenvalue problems. 
This modification of Calogero's method is especially useful for solving differential 
equations and eigenvalue problems with periodic solutions. Our motivation to use 
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trigonometric polynomials to approximate solutions of differential equations is also 
based on the Second Weierstrass Theorem, which states that every continuous L- 
periodic function / : [0, L] — t- M can be approximated, with any desired accuracy 
in the supremum norm, by an L-periodic trigonometric polynomial pQ. 

To construct representations of differential operators, we use trigonometric in- 
terpolation. Let u : [0, L) — 7- M be a continuous function. Partition the interval 
[0, L) into 2n + l points {xq, Xi, ...X2n\ such that = xq < xi < ... < X2n < L. Let 
Uj = u{xj) for all j G {0, 1, 2r2}. We construct a polynomial T[u] G such 
that 

(2) T[u]{xj) = Uj for all < j < 2n. 
This polynomial is defined to be [H] 

2n 

(3) T[u\{x) = ^Ujtj{x), 

j=0 



where 




for all X G [0,L). 

The functions tj satisfy tj{xk) = Sjk, where 6 is the Kronecker delta. Moreover, 
each tj is a trigonometric polynomial of degree n, which follows from trigonometric 
product-to-sum formulas. 

We note that the set of functions {tjj^^Q is linearly independent. Indeed, sup- 
pose that X]^=o '^i^ii^) — for all x G [0, L) and some real numbers ao, ai, . . . , a2n- 
By substituting x = Xj into the last equation, we obtain aj = for all j G 
{0,1,..., 2n}. 
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Theorem 1. For a given set of points S = {{xo,yo) , {xi,yi) , {x2n,y2n)} G 
M^, such that xo,---,X2n ^ [0,-^) are distinct, there exists a unique L-periodic 
trigonometric polynomial f of degree n that passes through all the points in S. 

The existence of a trigonometric polynomial 

f{x) = c + ^ ^Ofc COS (^^^^ + sin (^^^^ ^ 

whose graph contains the set S follows from the previous discussion, and its 
uniqueness is a consequence of the Fundamental Theorem of Algebra applied 
to an algebraic polynomial /i/ : C — )■ C associated with /. More precisely, 
/ is recast as f{x) = Zlm=-n CmC^*™^, where i^ = -1, Cm = Ham - ibm), 
C-m = ^(flm + ibm) foT m G {1,2, ...,n}, and cq = c. The polynomial hj is 
defined by h{z) = 2" YlZi=-n^rnZ'^ ^ where the complex variable z = e^*^. If g is 
another trigonometric polynomial that passes through the points in then the 
algebraic polynomial Q = hf — kg of degree at most 2n has 2n + 1 distinct zeros 

27r ■ 

Zk = e^*^*=, k = 0,1, . . . ,2n. By the Fundamental Theorem of Algebra, Q = 0, 
which implies f = g- 

Corollary 1. Let f be an L-periodic trigonometric polynomial of degree n, and 
let the projection T[f] be constructed using a partition {xj}^"g of [0,L), according 
to formula Then T[f] = f. Moreover, the system of functions {tj}^^^ forms 
a basis of ■ 

It is known that the set of functions 

/27r \ /27r \ f '^t^ \ 27r 

1, cos —X , sm —X , . . . , cos n—x , sm n—x 



L J \L J \ L J' \ L 

forms a basis of ■ We denote this system by {/i}i"o' where /o(a;) = 1, = 
cos (^x) , f2{x) = sin (^x) , . . . , f2n{x) = sin (^^nx) . We will use both bases in 
our discussion. 
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3. Matrix representations of differential operators via 
trigonometric interpolation 

In this section we construct finite-dimensional representations of linear differ- 
ential operators. We denote the space of all s times continuously differentiable 
functions on [0,L] by C^^j. 

Theorem 2. Let A be a linear differential operator of order s with a domain 
D{A) D C'pL]- -^^^ ^ ^ -^('^) {3^0,3^15 ■■■iX2n} he a partition of [0,L). Let 

u = {uq,ui, ...,U2nY , where Uj = u{xj) for all j G {0, 1, 2n}. There exists a 
(2n -|- 1) -dimensional square matrix A such that 

2n 

(5) T [AT[u]] (x) = J2 [^^]Mx), X G [0, L). 

j=0 

Moreover, 

(6) Ajk = Atk{xj) 
for allj,k G {0, 1, . . . , 2n}. 

Proof. To find a formula for the components of A, we substitute u = t^ into eq. (|5]), 
where A; G {0, 1, ... , 2n}. From Corollary [l| T[tk] = tt and eq. ([s]) becomes 

2n 

T [^tfc] (x) = ^ [Aefc] . tj (x) , X G [0, L) , 
i=o 

where {e^j^^Q is the standard basis of M^""'"-'^. 
By definition ([3]), 

2n 2n 
^ [^tA:(Xj)] tj(x) = Ajktj{x), X G [0, L). 

3=0 j=0 

Therefore, Atk{xj) = Ajk for all j, k G {0, 1, ... , 2n}. Let us verify that eq. ([s]) 
holds for an arbitrary function u G D{A). With u defined as in the hypothesis, 
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we obtain 



T [AT[u]] (x) =T 



2n 



A ^ Uktk 



2n 



k=0 
2n 



(x) 



j=0 lk=0 
2n T 2n 

j=Q \_k=0 
2n 

■■Y,[A%t,{x), xe [0,L), 
i=o 



as required. 



□ 



Definition 2. If a matrix A G ]R(2"+i)x(2n+i) g^j^^ linear differential operator 
A with domain D(A) D ^[o l] related by eq. ([s]), we say that the matrix A 
represents the differential operator A in ]R(^"+^). 

Let D he a {2n + 1) -dimensional square matrix that represents the operator ^ 
and let X be a {2n + l)-dimensional square matrix that represents the operator of 
multiplication by x. Using theorem |2| let us find explicit formulas for the entries 
of matrices D and X. For convenience, we define 



(7) 



2n 

n 

fc=0 
k^m 



sm 



^{Xm - Xk) 

L 



for all m G {0, 1, . . . , 2n}. For every vector v G M^, Diag(f) denotes the A^- 
dimensional square matrix whose diagonal is f; In denotes the A^-dimensional 
identity matrix. 
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Theorem 3. Let D,X be {2n + 1)- dimensional square matrices that represent the 
operator -p and the operator x of multiplication by x, respectively. Then 



D. 



IT 



L sin ^i^=L^ 
)^cot 



L 



ifj ^ m, 
zfj = m 



k=0 



and X = Diag(xo,a;i, . . . ,X2n)- Moreover, if P{x) = ao + aix + . . . + a^x* is an 
algebraic polynomial of degree s with real coefficients {aj}|=0) then the operator 
P (^) is represented by the matrix P (D). 

The proof of this theorem is a straightforward apphcation of formulas ([s]) and 
We note that if (7 is a trigonometric polynomial of degree n, then ^ is also a 
trigonometric polynomial of degree n. The representation of the identity operator 
1 is easily verified to be hn+i- 



4. Rank of matrix representations of differential operators 

Let P{x) = ao + Oi\x + . . . + a^x* be a polynomial of degree s and let ^ = P(^) 
be a differential operator with domain D[^A) D C'joL]- By theorem jsj the matrix 
P{D) represents the differential operator A on 1^^"+^. The main idea of our method 
is to approximate differential equation Au = f, where / G C*[o,l]5 by the system 
of linear equations P{D)u = f, where f = (/(xq), /(xi), . . . ,/(x2„))"^. To apply 
this idea to solution of boundary and initial value problems, we need to study the 
rank of the matrix P{D). The following theorem gives a complete answer to this 
question. 

Theorem 4. Let P be a polynomial of degree s given by 



(9) 



P{x) = ao + (^ix + . . . + cxsX^ , 
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where oq, ai, . . . ,as G M and Og ^ 0. Define the polynomials 

(p{k) = ao - a2 \—k j +a4\—kj - . . . + {-1} aq \—k 



and 



where g = 2 [|J is the largest even integer that does not exceed s, r = 2s— g— 1 is the 
largest odd integer that does not exceed s, 2i = q mod 4, and 2m = (r — 1) mod 4. 
Let D be the {2n + 1) -dimensional square matrix, defined by (|8|, that represents 



the operator ^ in M?""^^. Suppose that the syst 



em 



(10) 



m = 0, 

Cik) = 

has exactly m > solutions that belong to the set {1, 2, . . . , ra}. Then 



(11) rankP(D) = 2n + \ signaol — 2m 

Note. The results stated in this theorem do not depend on the choice of partition 
{xj}f"o of [0) L) with which the matrix D is constructed. 

Proof. We prove the theorem for the case ao 7^ 0. The case = can be proven 
in a similar fashion. 

By the rank-nullity theorem, it suffices to show that the dimension of the kernel 
of P{D) is 2m. Define the differential operator A = P{-^) and let u e ]R^"+^ We 
observe that P[D)u = if and only if Yl'i=o ~ 0' which happens 

if and only if AY^'^lo^i^ii^) ~ 0? where x G [0,L). Therefore, we will study the 
kernel of the differential operator A with its domain restricted to 7^^. 

Let g G be given by 

g{x) = b + ^Cfc cos ^^^^^ + '^k sin ^^^^^ ^ ' x G [0, L). 
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We will prove that Ag = if and only if 6 = and Ck = dk = for all k G 
{1, 2, . . . , n} that do not solve system ( [Io| ). That is, ker ^It^l = {0} if m = and 
keiAl-j-L = |cos (^fcjx) , sin (^fc^x) }^_^ if m > 1 and ki, A;2, • • • , /Cm are all the 
solutions of system ( [Io| ) in {1,2,..., n}. 
Indeed, because 



Ag{x) =bao + ^ 

k=l 

=bao + ^ i 

k=l L ^ 



^ d' /27r \ ^ ^ d' . /27r, 
Ck oii—^ cos I — fcx J + dk 2_^ai—^sm\^ —kx 



i=0 



i=0 



0(fi;)Cfc + -j-dkQ[k) I cos I -j^^^ 



(12) + 
= implies 6 = 0, 

(13) 



■^—CkC[k) + dk(p[k) I sm I -^"^a; 



27r/c 

Cfc0(fc) + —j-dkCik) = 0, and 



X G [0,L), 



(14) 



27rfc 



CkC{k) + dk(t){k) = 



for all A; G {1,2,..., n}. Suppose that A; G {1, 2, . . . , n} does not solve system (10) 
Then 0(/c) and ({k) are not both equal to zero, and 



(15) 



det 



0(A:) ^C{k) \ f27ik 



2iTk 



cik) m 



C\k) > 0. 



Therefore, the linear system ( I3p4 ) has the unique solution Ck = dk = 0. 

On the other hand, if is a linear combination of {cos (^^kjx) , sin {^kjX^ }™_^ 



or (yf = 0, then Ag = by formula (12). 



By a previous observation, the vector u G ker P{D) if and only if A Yli2o ^i^ii^ 



for all X G [0,L). If m = 0, then system (10) has no solutions in {1,2, . . . ,n} 



which implies YlfZo'^iti = and u = 0. In this case kerP(D) = {0} and 
rankP(D) is full. 
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If m > 1 and the numbers ki,k2, ■ ■ ■ ,km are all the solutions of system ( [Io| in 
{1,2,..., n}, then Y^f2o ^iti{x) is a linear combination of {cos (^'^kjx) , sin i^^kjx) }™_^- 
This occurs if and only if u is a linear combination of the vectors {vj, Wj}^-,^ C 
]^2n+i (jggned by the equations 

2n 



COS [^kjx \ = X][vj].ti(x) 



i=0 

and 

2n 



sin ij-kjx \ = ^ [wj].ti(x) 

^ ^ i=0 



for all j G {1,2,..., m}. The system of vectors {vj, Wjjjl-^ is linearly independent 
because the system of functions {cos (^/cjx) ,sin (x^i-^) }j=i linearly indepen- 
dent. Therefore, kerP(Z}) = span {vi, wi, V2, W2, v^, w^} and dim(kerP(D)) 
2m. 

□ 

Corollary 2. Consider the differential operator Q = (^)^ ^ + rn'^ where m G 
{1, 2, . . . , n}. The representation {^Y D"^ + rri^I of Q in ]R^"+^ has 

rank | ( + m^J | = 2n - 1, 



where D is the (2n + 1)- dimensional matrix that represents ^ in ]R^"+-'^ and / = 

Proof. Using theorem [i] with P{x) = + '^^j we compute (j){k) = — k"^ 

and C(^) = 0. System ( [Io| has exactly one solution, /c = m, in the set {1,2,..., n). 
Thus, rank ((^)^ + m^/j = 2n - 1. □ 



Corollary 3. Let D he a {2n + \) -dimensional square matrix representing ^ 
M2«+i. r/ien rank(DP) = 2n for all p G N. 

Proof. The result follows from the application of theorem |4] to the polynomial 
P{x) = x^. Note that ao = 0, one among {(f){k), ((k)} equals ± {^ky, and the 
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other equals zero. Because system (10) has no solutions in the set {1,2, . . . ,n}. 



rank P{D) = 2n. □ 

5. Numerical tests 
5.1. Boundary value problems. 

Example 5.1. The aim of this example is to test accuracy and convergence rate 
of the trigonometric interpolation method. Consider the boundary value problem 



(16) 



u" + u = 0, 

m(0) = u{tt/2) = 1, 



for a function u : [O, |] — t- M. Let a = {xi}^2Q be a partition of [0, 27r) such 
that Xfc = f for some A; G {0, 1, . . . , 2n} and let D,I e M(2«+i)x(2n+i) ^j^g g^j^^. 



dimensional representations of operators and 1, constructed using this partition 
and the period L = 271. We approximate the differential operator Q = ^ + 1 with 



domain C^[0, 2tx\ by + / and reduce problem (16) to the problem of finding a 
vector u G M^"+^ that satisfies 



(D2 + /)u = 0, 
Uo = Ufc = 1. 



(17) 

By Corollary [2| the matrix + / has rank 2n — 1, so the space of solutions 



of the first equation in (17) is two-dimensional, which implies the uniqueness of 



the solution of system (17). To solve system (17) in MATLAB, we define a new 
{2n + 3) X (2n + 1) matrix C whose upper block is + / and lower block consists 
of two rows ei,efc G 1^^"+^ where {eij'^Il'l^ is the standard basis of ]R^"'+-'^. Let 
f G M2'^+3 be defined as f = (0, 0, . . . , 0, 1, 1)' . We are now looking for the solution 
u G ]R^"+^ to a system of linear equations Cu = f. We used MATLAB to solve 
this system and compared the solution u with the vector of values of the exact 



solution u{x) = cosx + sinx of (16) at the = 2n + 1 points of the partition a. 
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The values of the error 



^max(iV) 



max 

ie{0,l,...,2n} 



Uj U [Xj^ 



are given in table [T] 

Table 1. The errors E^i^^{N) obtained using trigonometric inter- 
polation, Lagrangian interpolation, and the shooting method for 



boundary value problem (16), where N is the number of equidis- 



tant partition points of [0, 2tt). 



Method 


N 


5 


7 


9 


Trigonometric Interp. 




6.9389e-16 


3.7748e-15 


3.9968e-15 


Lagrangian Interp. 


-^max(-^) 


1.1466e-004 


5.2463e-006 


2.5575e-007 


Shooting Method 




7.3e-3 


3.2e-3 


1.8e-3 



We compared the accuracy of the trigonometric interpolation method described 
in this paper with the accuracy of the Lagrange interpolation method [21 15] and the 
standard shooting method. To apply the shooting technique, we recast boundary 



value problem (16) into two initial value problems: 



V" = -V, 

v{0) = l,v'{0) = 



and 



w = —w, 
w{0) = 0,w'{0) 



1, 



which we solve by the method of finite differences. The solution to problem (16) 
can be found by the formula 



u 



V + 



l-v( 



-w. 



Trigonometric interpolation is a suitable choice for problem (16), since its solution 
is 27r-periodic. In fact, u{xi) = Ui for all z G {0, 1, . . . , n} if u is the exact solution 
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of Cu = f and the errors in the first row of table [l] are computational errors that 
occur due to the limited number (sixteen) of significant digits of each component 
of u stored in the computer memory. Table [l] shows that in this case our method 
yields significantly higher accuracy when compared to the shooting method or a 
similar method that uses Lagrangian interpolation [21 [5] . 

Example 5.2. Consider the boundary value problem 

u" + u' + u = 0, 
u{0) = u{^^) = 1. 

Let D be the finite representation of ^ given by formula ([s]), where L = By 
Corollary |2| the rank of + D + / is full, so we make the substitution 

u{x) = X — v{x) + 1 



to eliminate the boundary conditions in (18). The substitution reduces prob- 



lem (18) to the differential equation 



(19) [(2 -a){x + l) + x^]v + (-2a + Ax - ax + x^) v + {-ax + x^) v" = -1, 
where a = with the conditions that require f (0) and v{a) to be finite. 



To approximate problem (19) with a system of linear equations using trigono- 
metric interpolation, we must choose the period L of the trigonometric functions 
used in the approximation scheme. Although a natural choice seems to be L = 
it would impose an implicit condition v{0) = v ^^^^ Having this in mind, we 
double the length of approximation interval and use the period L = in for- 
mula ([S]). We partition the interval 0, , into {2n + 1) points Xj = ■ ^, 
where j G {0, 1, ... , 2n}, and consider this as a partition of [0, -^). We compute 
the matrix D using formula ([S]). To improve the conditioning of a matrix repre- 
sentation of -|- ^ "'"-'-)' matrix D = \['~^Z)\E' instead of where 
^ = Diag{ipo, ^i, . . . , ip2n) (read Theorem [g] for more insight on the use of D). 
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We tlien solve the system of linear equations 

[(2 - a){X + I) + + (-2a/ + AX - aX + X^) D 
+ {-aX + X^) D'^]v = - diag / 

for V G ]R^"+-^, and calculate the vector u G M^""'''^ using formula Uj = Xj — j '?AiVj+ 
1, where i G {0, 1, . . . , 2n}. The vector u approximates the values of the exact so- 
lution 



u[x) = exp 



V2 



exp 



4x/3 



1 I sm h cos 



to problem (18) at the partition points. 



In the following table we summarize the approximation error 

-Emax(iV) = max \Ui-u{Xi)\ 

j6{0,l,...,2n.} 

for several values of = 2?2 + 1 . 

Table 2. Values of the error -Emax obtained when solving prob- 



lem ( 18 ) using trigonometric interpolation. 



Number of Partition Points N 


5 


11 


21 


^max(iV) 


5.6458e-4 


3.7868e-6 


3.7706e-9 



Example 5.3. Consider the boundary value problem 

(^)^(^) + 13(^)^(^) + 36n = 0, 



(20) 



m(0) = u{l) = 1,m'(0) = u'{l) = 0. 
We partition the interval [0, 1) by points {xj}^"q and approximate the differential 



equation in (20) by the system of linear equations Au = for u G M , where 

2 

/ I \ . / I \ 

A 



and D is computed by formula ([8| using the period L 
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By theorem |4| rank A = order A — 4. Therefore, we can complement the hn- 
ear system Au = with four additional equations that represent the boundary 



conditions in (|20|). More precisely, we solve the linear system 

Ail 



0, 

UO = U2n = 1, 







or, equivalently. 



where the block 



B 



( 1 













\--\ 


;:) 





















0,0 



^0,1 



a 



0,2n-l 



\ 

1 

Do,2n 



\ D2nfl D2n,l ■ ■ ■ D2n,2n-1 -D2n,2n / 



and b = (1,1,0,0)^. 



We compute the approximate solution of (20) by the formula 

2n 

(21) Mapprox(a;) = ^UA(X). 



i=0 



In fig. [T| we graph the approximate solution Mapprox obtained using the partition 



X'l 



2n- 



-j, where z G {0, 1, . . . , 2n} and n = 5. 



5.2. Eigenvalue problems. Eigenvalue problems play an important role in Physics 
and Applied Mathematics. Solutions to many partial differential equations that 
model heat conduction, wave propagation, and other physical phenomena, can be 
reduced to the solution of a Sturm-Liouville eigenvalue problem. 

In this section, we present solutions to several Sturm-Liouville eigenvalue prob- 
lems obtained using trigonometric interpolation. Below we state some background 
information useful for understanding applications. 
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Figure 1. The approximate solution to (20) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Theorem 5. Let A = ao+ai ^ + • • -+0^5^ be a differential operator with domain 
D{A) D C^Q^, where ao, ai, . . . , G M and 7^ 0. Let (a, /) G M x D{A) be an 
eigenpair of the operator A. If f , then the finite- dimensional representation 
A = af) + aiD + . . . + a^-D* of operator A, constructed using a partition {xjj^^Q of 
[0,L), has an eigenpair (q;,w) G M x M^"'+-'^ such that Yld=o^i^ii.'^) — /(^) /^'^ '^^^ 
X G [0,L). 



This theorem was proven by Calogero [7] . Because we use different notation, we 
supply a short proof in appendix |X] for the reader's convenience. 

Let \Ef = Diag('?/'o; V'l) • • • i'4'2n), where ipj^s are defined by formula (jT]). In the 
following examples, we use the matrix D = "^^^D"^ instead of D to solve eigenvalue 
problems because it yields representations of differential operators that are usually 
better conditioned. 
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Theorem 6. // {a, /) G M x is an eigenpair of a differential operator A defined 
as in theorem \^ then the matrix A = aQ + ait) + . . . + ast)^ has an eigenpair 
(a, w) G M X M2n+i^ y;/ie^e E?=o ^^^^^^(a;) = f{x) for all x G [0,L). 

This theorem is a direct consequence of Theorem [s] and the formula D = \l/£)\l/~^. 

Example 5.4. This example aims to test the accuracy and convergence rate of 
the trigonometric interpolation method applied to eigenvalue problems. Consider 
the eigenvalue problem 

' -y" = \y, 

y(0)=y(l)=0. 

Let {xj}^"g be a partition of [0,1). Let the matrix D = \E'~^D\1', where D is 
defined by eq. (j8| with L = 1. Let X = Diag (xq, . . . , X2n) and / = hn+i- To get 



(22) 



rid of boundary conditions in (22), we make a substitution y{x) = x{x — l)v{x). 



Via this substitution, the eigenvalue problem (22) is recast as 



- [(x^ - x)v" + {Ax - 2)v' + 2v] = Xx{x - l)v, 

with the condition that requires v{0) and f (1) to be finite. The last eigenvalue 
problem is then approximated by the matrix eigenvalue problem 



(23) -(^(X2-X)i)2 + (4X-2/)Z) + 2/)v = A(X(X-/))v, vG 



[)2n+l 



We use MATLAB to find the eigenvalues of (23) using the partition Xj = ^ for 
i G {0, 1, . . . , n — 1}. The exact eigenvalues for this system are known to be vr^j^ 
for j G N, and we compare the approximate eigenvalues with exact eigenvalues 
in table [3} which also illustrates the convergence of this scheme as the number of 
equidistant partition points increases. 



Using the trigonometric interpolation approach, we solved several eigenvalue 
problems and compared the approximate eigenvalues we obtained with the lower 
and upper bounds found in [9]. 
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Table 3. Approximate eigenvalues of problem (22), where Aj de- 
notes the 2*'^ eigenvalue. 



Number of partition points 


Ai 


As 


All 


Al5 


19 


9.906 


247.699 


1200.521 


2240.381 


31 


*.883 


**7.088 


*196.079 


**24.557 


51 


*.*75 


**6.867 


***4.858 


***1.891 


101 


* **-|^ 


*** 772 


**** 2g 


***0.956 


Exact eigenvalue 


9.8696 


246.740 


1194.222 


2220.660 



Example 5.5. Consider a genetics eigenvalue problem 
(24) 



1(1 - x^) {-u" + 2a [{xq - x) u]') = Aw, 
m(-1) = m(1) = 0, 



where xq, a G M. 



Using the substitution u{x) = (1 — x'^)v{x), we approximate (24) by the matrix 
eigenvalue problem A\ = Av, where v G M^""*"^ and 



(25) 



4 ^ ^ 



X + -a{l- X^) {xol -X) b 



+ ^ + aX (X - xd) -\a{l- X^) 



Results in tables |4] and |5] are obtained using the partition Xi 
i G {1, 2, . . . , n}, where n = 25 and the period L = 4. 



2i-n-l 



for all 



Example 5.6. Consider the eigenvalue problem 



(26) 



u" + A (1 + sin x) ti = 0, 
m(0) = M(7r) = 0. 
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Table 4. Aj denotes the i^^ eigenvalue of problem p4|), where a 



and Xo = 0.6. Upper and lower bounds for Aj are taken from [9]. 





Approximate Eigenvalues 


Lower Bound 


Upper Bound 


Ai 


0.3938 


0.3931 


0.3939 


A4 


4.9514 


4.85 


4.96 


As 


17.9515 


16.7 


17.96 


Aio 


27.4552 


24.6 


27.46 



Table 5. Aj denotes the r eigenvalue of problem (24), where a 



2, and Xq = 0.4. Upper and lower bounds for Aj are taken from [U]. 





Approximate Eigenvalues 


Lower Bound 


Upper Bound 


Ai 


0.2991 


0.2982 


0.2992 


A4 


4.9625 


4.73 


4.97 


As 


17.9570 


15.3 


17.96 


Aio 


27.4585 


21.7 


27.46 



We make the substitution u{x) = y{l — y)v{y), where y = x/n. We approximate 
problem (26) by the matrix eigenvalue problem = Av, where v G IR^"'^^ and 



A = ^{I + sin(7rX))~^ (-D^ + 2[X{I- X)]'^ \{2X - I) D + 1 



The matrices X and D are constructed using a partition {xj}|"Q of [0,2). We 
use the interval [0, 2) instead of [0, 1) for approximation of the function v to avoid 
imposing the implicit condition v{0) = v{l). 

for j G {1,2, . . . ,n}. 



Results in table [6| are obtained using the partition Xj 
n = 21, and the period L = 2. 



n+l 
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Table 6. Aj denotes the i^^ eigenvalue of problem 26 Upper and 
lower bounds for Aj are taken from [9]. 





Approximate Eigenvalues 


Lower Bound 


Upper Bound 


Ai 


0.5403 


0.540282 


0.540319 


As 


5.4486 


5.411 


5.449 


As 


15.3126 


14.54 


15.313 


Ae 


22.0972 


21.39 


22.1 
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Appendix A. Proof of theorem [5] 

Proof. Let (a, /) G M x be an eigenpair of A. Choose a partition {xjjf^Q of 
[0, L) and let A = + aiD + . . . + UgD^, where D is given by formula ([s]). Let 

f = ifiXo), fiXi), . . . , f{x2n)f- 

Because / G T^, Af G Tn[x] and the equation Af{x) = af{x) is equivalent to 
the equation 

2n 

by corollary hi The last equation is equivalent to Af = af. Hence, (a,f) G 



2n 



j=0 



X) 



XG [0,L), 



t)2n+l • 



is an eigenpair of A. Moreover, f{x) = Yll^o^i^ii^) ^ ^ [0'-^) 



required. 



□ 
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